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Following the model independent approach to deuteron photodisintcgration with linearly 
polarized 7— rays, we show that the measurements of the tensor analyzing powers on 
aligned deuterons along with the differential cross section involve five different linear 
combinations of the isovector El J v ;j = 0,1,2 amplitudes interfering with the isoscalar 
Ml„ a nd E 2 a amplitudes. This is of current interest in view of the recent experimental 
finding E21 that the three E1 J V amplitudes are distinct and also the reported experimental 
observationSSlon the front-back (polar angle) asymmetry in the differential cross section. 



1. Introduction 

The earliest studies on deuteron photodisintegration and its inverse reaction viz., 
n-p fusion go back to more than seven decades ^4 These reactions are relavent 
to nucleosynthesis scenarios ^ of cosmological interest from the Big Bang to stel- 
lar evolution under various environments. The earliest estimates of the reaction 
rates by Fowler, Caughlan and Zimmermann ^ used the theoretical calculations 
^ of deuteron photodisintegration normalized to the then available thermal neu- 
tron radiative capture cross section measurements Since the deutron is loosely 
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G.V.K. Academy, Bangalore, India 
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bound, it is destroyed in stellar environments even before the stars reach the main 
sequence and as such it has been observed EI that 'the ratio of the primordial abun- 
dance of deuterium to that observed today could be anywhere between 1 and 50'. 
Although the measurements ^ of deuterium abundance in high-red-shift hydrogen 
clouds have improved the scenario, Buries et al 121 have highlighted the need for the 
precise knowledge of the reaction rates for d + -f^n + p&t astrophysically relevant 
energies which correspond to the c. m. neutron energy range 25 to 200 keV. 

While studies on photodisintegration of the deuteron are well documented 
at lab photon energies > 2.33 MeV and a large number of measurements exist 
for the n-p fusi on cro ss section at thermal neutron energies which correspond to 
E% = 10~ 5 keV^^ES there are only two measurements at astrophysically relevant 
energies by Suzuki et al.,^21at lab neutron energies E% = 20, 40 and 60 keV and by 
Nagai et al.^at E% = 550 keV as shown in Table 1. It was known quite early that 
the thermal neutron cross section is dominate d b y the isovector magnetic dipole 
amplitude denoted as Ml v . Breit and Rustgi ^ were the first to propose a po- 
larized target-beam experiment to look for an isoscalar magnetic dipole amplitude 
denoted as Ml s , in view of the then existing 10% discrepancy between theory and 
experiment. Althoug h th e meson exchange current contributions (MEC) suggested 
by Riska and Brown ^3 explained this discrepancy in the cross section with sur- 
prising accuracy, the measured values for analyzing powers in p(n,j)d as well as 
f° r n eutron po larization in photodisintegration of the deuteron were both found to 
differ I I ^ 1 u l from theoretical calculations which included MEC effects. In fact, it 
has also been pointed out by Rustgi, Vyas and Chopra ^] that there is an unam- 
biguous disagreement between theory and experiment when a comparison is made 
with the data at a photon energy of 2.75 MeV. They pointed out moreover, that 
the inclusion of the two body effects widens this disagreement. 

Several theoretical studies have been reported on d + 7 ^ n + p based on 

a variety of deuteron models, including those at quark level. Some of these cal- 
culations take into account meson exchange and other currents. Nagai et al., ^Sl 
have quoted the theoretical results of Sato et al.,SH which included MEC, isobar 
currents and pair currents. These theoretical results show that the Ml v strength 
decreases sharply with increasing energy, while El v i.e., the isovector electric dipole 
amplitude starts appearing in the energy range of interest to astrophysics and in- 
creases with energy there after. Moreover, the theoretical calculations quoted by 
Nagai et. al., show that the El v and Ml v strengths are approximately equal 
at around cm neutron energy E n of order 500 keV. The relative strengths of ML, 
and El v have also been discussed in several studies using effective field theory I^SI 
While the role of isoscalar Ml s have been discussed by Chen et al.,^Sl an d Park 
et al., 123 using different versions of effective field theory. The predictions of these 
two ver sions led to an experimental measurement of the 7 anisotropy by Muller et 
al., ^S! ; w hich however could not decide in favor of either of the approaches. The 
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calculations of E2 S contributions have also been discussed by Chen et al.,^land 
Hadjimichael et al.,^l Although the El v transition dominates photodisintegration 
of deuteron, the E2 S contributions, though small, has been considered to be still 
significant "01. 

An observable which is sensitive to the presence of isoscalar Ml s and E2 S transi- 
tions from the triplet S-state is the circular polarizatio n of the emitted radiation 
with initially polarized neutrons. The first measurement to detect the presence 
of isoscalar amplitudes was not quite encouraging but a subse quen t measurement 
&i yielded a value P 1 = -(2.29 ± 0.9) x 10~ 3 . An attempt 33 to explain the 
large measured value by introducing a six quark admixture in the deuteron wave 
function led, however, t o a disagreement with the precise deuteron magnetic mo- 
ment. Later calculations ^21 in the zero range approximation and the wavefunction 
for a Reid soft core potential led to a theoretical prediction P-y of the order of 
—1.1 x 10~ 3 with an estimated accuracy of 25%. The latest measured value of 
Ry = — (1.5±0.3) x 10~ 3 is in reasonable agreement with the theoretical calculation 
The importance of measuring the photon polarization with initially polarized 
neutrons incident on a polarized proton target has been pointed out When the 
initial preparation of the neutron and proton polarizations P(n) and P(p) are such 
that they are either opposite to each other or orthogonal to each other, the inter- 
ference of the small isoscalar amplitudes with the large isovector amplitude could 
substantially contribute to the observable photon polarization. Attention has also 
been focussed on the spin response of the deuteron, i.e., the asymmetry of photoab- 
sorption with respect to parallel and antipa rallel s pins of photon and deu teron and 
the Gerasimov-Drell-Hearn (GDH) sum 

rule [37[38| by 

Arenhovel et al. , who have 
also considered other polarization observables in electromagnetic deuteron break- 
up SOI anc j compare their formalism with that of Dmitrasinovic and Gross SB. We 
might also mention an experimental measurement of recoil proton polarization in 
electrodisintegration of deuteron by polarized electrons S21 and a theoretical discus- 
sion of photodisintegration of pol arize d deuterons by unpolarized photons S21 using 
the model independent formalism ^1 

As the Big Bang Nucleosynthesis (BBN) entered the precision era^Sl, an accurate 
estimation of the primordial abundance of deuterium became important and it has 
been referred to as the Cosmic Baryometer. The work of Buries et al., ^ catalyzed 
the first experimental study SSI on the analyzing power 2(0) in deuteron photodis- 
integration at El = 3.58 MeV and 9 = 150°, using 100 % linearly polarized pho- 
tons from the High Intensity Gamma-ray source (HI GS) at the Duke F ree Electron 
Laser Laboratory. It was followed by further studies B7 | 48 | 4 9 50 ol 52] These mea- 
surements have been analyzed, using a theoretical formalism^! where the Ml v and 
El v multipole contributions (which are dominant at energies of thermal neutron 
capture and deuteron photodisintegration res pect ively) were calculated separately. 
The model independent theoretical formalism f or deuteron photodisintegration 
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Table 1. Measured cross-sections 





for p(n, j)d 






for d(7, n)p 




Energy 


Cross-section 


Reference 


Energy 


Cross-section 


Reference 


E% keV 


nib 




E% MeV 


mb 




10" 5 


334.2 ±0.5 


EH 


2.33 


0.683 ±0.053 ±0.042 


EH 


10~ 5 


332.6 ±0.7 


M 


2.52 


0.983 ±0.039 ±0.061 


E2 


20 


(318 ±2) x 10~ 3 


M 


2.62 


1.300 ±0.029 


El 


40 


(203 ± 19) x 10~ 3 


M 


2.754 


1.456 ± 0.045 


an 


G4 


(151 ± 7) x 10~ 3 


M 


2.76 


1.474 ±0.032 


El 


550 


(35.2 ±2.4) x 10" 3 




2.79 


1.47 ±0.03 ±0.09 


M 



with linearly polarized photons showed that the differential cross section in d(ff, n)p 
contains a term representing the interference between the El v and Ml s amplitudes, 
which is non zero if the three El v am plitudes EV V , j — 0, 1,2 are unequal. In a 
recent publication, Blackston et al., EShave reported the first experimental obser- 
vation of the splitting of the three El v amplitudes. 

Therefore, the purpose of the present paper is to employ the model independent 
theoretical approach EH to carry further the study of the contributions of isoscalar 
Ml s and E2 S amplitudes in photodisintegration of deuterons by linearly polarized 
photons using aligned deuteron targets. 



2. Aligned deuteron 

The deuteron is a spin-1 nucleus, which gets polarized in different ways in different 
external environments. For example, an oriented deuteron target can be pro- 
duced under the influence of an external uniform magnetic field, whose direction 
determines the axis of orientation. On the other hand, when the deuteron is exposed 
to external electric quadrupole fields generated by surrounding electrons in crystal 
lattice sites ^1, the spin of the deuteron is aligned. States of polarization of spin-1 
nuclei exist, which are more complex and are in fact multiaxial ^1 The state of 
polarization of a spin-1 nucleus is, in general, specified by the spin density matrix 

rp (1+1*5 + £3 |(*-a+*-l) V3* 2 _ 2 ) 

-f*i+ £ ? 1-^*0 §(£i+£i) (1) 

[ v/3*i -m+ti) i-fij+^j 

which is hermitian and is parametrised in terms of the Fano statistical tensors 
tg,q — k,k — l,...,—k with k = 1 referring to its vector polarization and k — 2 
referring to its tensor polarization. The normalization used for the in eq. (1) 
follows I 56 | 58 | 59 [60| rows an( ^ co i umns Q f the above matrix are labeled by the 
states \lm > of the deuteron spin-1 with magnetic quantum numbers m = ±1, 0, —1 



i 



i 



March 6, 2013 18:49 WSPC/INSTRUCTION FILE 

aligned ' deuteron' ij mpe 



Photodisintegration of aligned deuterons at astrophysical energies using linearly polarized photons 5 

in that order w.r.t any chosen quantization axis, Q referred to usually as the Z-axis. 
The target is said to be unpolarized if all the = 0. Defining vector polarization 
through 

*o = \J\ p z t±i = t\V*{P x + iPy) (2) 

the target is said to be purely vector polarized if the vector polarization F^O but 
tq = 0. Defining the tensor polarization through a traceless symmetric 2nd rank 
cartesian tensor T Q( g; a,j3 = x,y,z given by 

to = t±i = T— j=(T xz ±iTy Z ) t\ 2 = -^-j=(T xx — Tyy±2iT xy ), (3) 

a spin 1 nucleus like the deuteron is said to be aligned, if its vector polarization is 
zero and tensor polarization is non-zero. 

Although the energy levels of spin 1 nuclei like deuteron (with non-zero electric 
quadrupole moment) have been studied in external electric quadrupole fields 
Ramachandran, Ravishankar, Sandhya and Swarnamala Sirsi are probably the 
earliest to draw attention to the non-orientedness of the nuclear polarization in 
such an environment when the asymmetry parameter r\ of the external electric 
quadrupole field with strength A is non-zero. Precise estimates of the parameters 
A and r\ when they are embedded in various compounds and at different sites are 
available I 55 | 61 [ The estimates of i 2 are also available ^^31 a t temperature of order 
mK. With the present day technology, in the case of ultracold atoms in optical 



lattice, the lowest temperatures reached are of the order of nK and pK ES. R 
could therefore be feasible to prepare aligned deuteron targets with higher values 
of tensor polarization. 

An external electric quadrupole field te nsor V a p is represented by Vxx , Vyy and 
Vzz in the Principal Axis Frame (PAF) where V a @ cx 6 a /3. When a spin-1 nu- 
cleus with a non-zero nuclear electric quadrupole moment Q is exposed to such an 
external electric quadrupole field, the resulting polarization is describable in terms 
of the Fano statistical tensors t k q where = = Zm(t|) = 0, which defines the 
Principal Axes of Alignment Frame (PAAF) "21 which coincides with PAF in this 
case when no other external fields like magnetic field are present. The interaction 
Hamiltonian is well known and is of the form 

H mt = A[(3J 2 Z - J 2 ) + V (J 2 X - J Y )], (4) 

where Jx , Jy & n d Jz are the cartesian components of the nuclear spin J with 
J 2 = J ■ J and A is proportional to the nuclear electric quadrupole moment Q 
through 

A = -OV ZZl 7? = Vxx ~ Vyy ; \V zz \>\Vyy\>\Vxx\. (5) 
4 Vzz 
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Table 2. Polarization parameters of aligned deuterons at temperature T =0.1 mK 



Compound 


A 


-v 


t 2 


t 2 


acetamide (ND2) 


0.6242 


0.4022 


-0.4584 


-4.6791E-02 


chloroacetamide (ND2) 


0.6260 


0.3772 


-0.4601 


-4.3949E-02 


formamide (ND2) 


0.5752 


0.3712 


-0.4223 


-4.1353E-02 


phthalamide (ND 2 )ii 


0.6457 


0.3722 


-0.4749 


-4.4035E-02 


L-asparagine (ND2) 


0.6377 


0.3826 


-0.4688 


-4.4990E-02 


L-asparagine hydrate (ND2) 


0.6797 


0.3272 


-0.5010 


-3.9644E-02 



The resulting eigen states \X >, \Y >, \Z > with ^4(1 ± rj) and — 2A energy eigen 
values respectively are given by 

LY>=i=(|l-l>-|ll>); E X = A(1 + V ) 

\Y >=-J=(|l-l>+|ll>); E Y =A(1- V ) 

|Z>=|10>; E Z = -2A, (6) 

in terms of magnetic substates \lm > defined with respect to the Z-axis of the 
PAAF, where t 2 ±1 = 

The state of the aligned deuteron is characterized by two axes which are along 
n\ = (0i,(j>i) and ri2 = (^2,</ ) 2) in PAAF. The ±r?i and ±ri2 are determinable by 
setting t\ = 0, which leads to a quadratic equation having in general two solu- 
tions. In the particular case where ii\ and ti2 are collinear, the polarized spi n-1 
system is oriented. Otherwise it is said to be nonoriented. Some typical cases ^1 
are illustrated in Table 2. 



3. Model Independent Theoretical Formalism 

Following 1331 a nd using the same notations, the reaction matrix for d + 7 — > n + p 
with linearly polarized photons is 

1 s+l 

M = E E (S X (s,l)-F X (s)), (7) 

s=0 A=|s-l| 

where S^(s, 1) are irreducible tensor operators of rank A in hadron spin space 
connecting the initial spin 1 state of the deuteron with the final singlet and 
triplet states, s = 0, 1 of the n — p system in the continuum. The irreducible tensor 
amplitudes, F*(s) are given by 

^(s)= E -^M. (8) 

(U=±l 
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where the four irreducible tensor amplitudes, J-*(s,n) may be written, taking into 
consideration both the isoscalar amplitudes Ml s and E2 S , in addition to the isovec- 
tor Ml v and the three EV v ,j = 0, 1, 2 amplitudes, as 

jS(l l/ i) = ^l t (0)/ o °(l,l,rt; 
^(0,/x) = -iMl v ri(0,l,ti) 
^(l,M) = -^l K (l)/ I i(l,l, M ) + iMl s /i(0,l,/i); 
J*(l, /x) = ^1„(2)/2(1, 1, M ) + E2 s /2(0, 2, /i). 



0) 



The -E1„(A = 0,1,2) in the above equation are linear combinations of the EV V 
amplitudes with j — 0, 1, 2. They are explicitly given by 



(10) 



~Ei v {oy 




"13 5" 




~EV=°~ 


£1,(1) 




2 3-5 




El j=1 


_E1 V (2)_ 




.2-3 1 . 




EV= 2 

_ V 



and f„(l, L, /i) are given by 



/*(Z, L, n) = (2 7 r) 3 / 2 M 7r+ C(/, L, A; m u -/x, v)Y h 



(11) 



which take care of the angular dependence and also the dependence on photon cir- 
cular polarization [i = ±1. 



The differential cross section with linearly polarized photons is then given by 



where 



= ^—[a + 6sin 2 9(1 + cos 20) - ccosfll. 
ail 6 



c = 4V6[Re[(El v (l))Ml*} + Re[(El v (2))E2*]] . 



(12) 



(13) 



This result, viz., eq. (12), for the differential cross section with polarized photons 
bears comparison with the form for the differential cross section obtained by Rustgi 
et al.j in their approximations D and E which take into account the Ml s and E2 S 
amplitudes. The cos 6 term is present in both these approximations. The presence 
of the term with coefficient c can easily be detected by measuring the front-back 
(polar angle) asymmetry 



da ... da .„ 27r 2 
^W- ^0-0) = --ccos< 



(14) 



In the abstract of his experimental study of d(7, n)p at laboratory photon energies 
3.5 to 10 MeV, Sawatzky-^ has observed: "Analysis of these data revealed a striking 
discrepancy with the prevailing low energy models of the Nucleon-Nucleon (N-N) 
interaction. The most prominent feature of the disagreement lies in an observed 
front-back (polar angle) asy mmet ry in the cm. photoneutron yield that is not rep- 
resented in the theory ( e .g.lHMZ|) The magnitude of this discrepancy rises as the 
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7-ray energy falls towards the threshold" . The last sentence in the above quote is 
particularly relevant for the discussion of the problem at astrophysical energies. 

The asymmetry represented by eq. (14) could also be caused by /—waves, i.e., with 
I = 3 in the final state, but the electric multipoles in / waves do not interfere with 
the Ml s amplitude. Note also that c reduces to zero, if all the three El J v ;j = 0, 1, 2 
amplitudes are equal: if so, the asymmetry represented by eq. (14) could possibly 
be caused by an EV V =2 multipole amplitude leading to I = 3 (which is different 
from three El^,j = 0, 1,2 amplitudes leading to I = 1) interfering with the E2 S 
amplitude. But it was found recently ^ that that the three EV V amplitudes are 
not equal at 14 MeV and 16 MeV. Such an inequality is likely to be present even at 
low energies like 3.5 to 10 MeV employed in the experimental studies of Sawatzky 
SS. Therefore, it is worth while to carry out further studies on the inequality of the 
EV V amplitudes leading to p- waves at lower energies and their interference with the 
Ml s an d E2 S amplitudes. 

4. Analyzing powers for deuteron disintegration by linearly 
polarized photons 

Let us therefore consider the differential cross section for photodisintegration of 
aligned deuterons by linearly polarized photons, which is given by 

(^U gned = ^TrrMpMt], (15) 
where p is given by equation (1) with non-zero t 2 . This leads to 



(17) 



where the analyzing powers are given by 

A \ = ^ E (-!) A ( 2s + 1)[AP'M12 S A; l\')(F x (s) ® J** '(*))» 

which are expressible in the form, 

Al = V2tt 2 [oo + b sin 2 + 2c sin 2 cos 20 + d cos 0] (18) 



where 



a = [^\M1 V \ 2 - 2\M1 S \ 2 - 2\E2 S \ 2 - 3|£1„(1)| 2 - 3|S1,(2)| 2 

-l8Re[{El v (2)El* v {l)} - 6Re{E2*Ml*)] , (19) 

b - 2[3|£T„(1)| 2 + 4|£1„(2)| 2 + Re{2El%{Q) + 9El* v (l))El v {2)] , (20) 

c = [^|£1„(1)| 2 - ±\E1 V (2)\ 2 + 2Re(El v (2)El* v (0))} , (21) 
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d a = -l|fl e (£l„(l) + 3£71„(2))Mi: + ^e(3M(l) + El v (2))E2* s ]. (22) 



and 



and 



A\ = 2\/37r 2 [a 2 + b 2 sin 2 6 + c 2 sin 2 0e 2 ^ + d 2 cos 9] . (23) 

«2 = [||^1.| 2 - |M1 S | 2 - |£l s | 2 - \\EI V {1)\ 2 - ||M„(2)| 2 

+3Re(El v (2)El* v (l)) + Re(E2 s Ml*)] , (24) 

6 2 = [^\E1 V (1)\ 2 + \\El v {2)\ 2 - 2Re(El v (2)El* v (0))] , (25) 

ca = [ - |S1,(2)| 2 - i?e(2£i;(0) + 3£i;(l))£l„(2)] , (26) 



ft = ^Re[(El v (2) -M„(l))(Mi; -i^)]. (27) 

5. Discussion 

It is important to note that the coefficients of cos 9 in the differential cross section ^ 
as well as in the analyzing powers Aq and A 2 (i.e., c, do and ft) involve interference 
between the three isovector E1 V (X) amplitudes given by eq. (10) and the isoscalar 
Ml s and E2 S amplitudes. Using eq. (10) and expressing 

c = Re{aMl*+pE2* s ) (28) 

do = Re{a Ml* + faE2*) (29) 

ft = Re(a 2 Ml* + ft£2*), (30) 

where 

a = 4V6(2£LT° + 3EV V =1 - 5EV= 2 ) (31) 

p = AV&(2El{= - Si^ir 1 + EV= 2 ) (32) 

a Q = -4V6(E1{=° - 3-Bir 1 + 2£ir 2 ) (33) 

ft = -V6(Elj= + 3£lf° - 4£;ir 2 )- (34) 

«2 = -ft = yf (£ii=° - sir 2 ) (35) 

in terms of three EV v ,j = 0, 1, 2 amplitudes. 
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It is interesting to note that the above equations (31) to (35) involve five dif- 
ferent linear combinations of the three EV V amplitudes and as such facilitate their 
study individually at lower energies of interest to astrophysics. The experimental 
finding 1^2! at 14 MeV and 16 MeV that all the three EV V amplitudes are not equal 
is encouraging. Since the possible role of the M\ s and E2 S amplitudes have been 
discussed by several authors using different theoretical formalisms over the last 50 
years, we feel that it is desirable to carry out the measurements on the tensor ana- 
lyzing powers Aq and A\ in photodisintegration of aligned deuterons in addition to 
the differential cross section ^fe. Such experimental studies are expected to clarify 
purely empirically the role of the isoscalar Ml s and E2 S amplitudes at the range 
of energies of interest to astrophysics. 
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